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Abstract. In this paper we determine the constraints on cosmological parameters using the CMB data from the Wmap exper- 
iment together with the recent power spectrum measurement of the SDSS Luminous Red Galaxies (LRGs). Specifically, we 
focus on spatially flat, low matter density models with adiabatic Gaussian initial conditions. The spatial flatness is achieved with 
an additional quintessence component whose equation of state parameter w e g is taken to be independent of redshift. Throughout 
most of the paper we do not allow any massive neutrino contribution and also the influence of the gravitational waves on the 
CMB is taken to be n egligible. The analysis is carried out separately for two cases: (i) using the acoustic scale measurements 
as presented in Hiitsi (2006), (ii) using the full SDSS LRG power spectrum and its covariance matrix. We are able to obtain a 
very tight constraint on the Hubble constant: Ho = 70.8*?i km/s/Mpc, which helps in breaking several degeneracies between 
the parameters and allows us to determine the low redshift expansion law with much higher accuracy than available from the 
Wmap + HST data alone. The positive deceleration parameter q Q is found to be ruled out at 5.5cr confidence level. Finally, we 
extend our analysis by investigating the effects of relaxing the assumption of spatial flatness and also allow for a contribution 
from massive neutrinos. 
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1. Introduction 



Since the flight of the Cobe [j satellite in the beginning 
of 90's the field of observational cosmology has witnessed 
an extremely rapid development. The data from various 
C osmic Microwave B ackgrou nd (CMB) experim ents (Wmap 
| (Benn ett et all l2003b . Cobe dSmpot et al J 1 1 992h. Archeop s 
3 (Beno itetalJl2003l). Boom erang Fl ( Netterfi ekTet al.ll2002h . 
MAxiMApldHananv et aljiooob. CBiRdPearsqn et alj2003l) . Vsa 



Q dScott et alj2003h . DasiFI ( Halverso net alj|2002l) etclsuper- 
no ya surveys (ScpB ( Perlmutter et al. 1999), High-Z SN Search 
( Riesset all 1998b ) and large galaxy redshift surveys (SDSS 



dYork et all2000l) . 2dFGRS[ n l fcolless et alfeOOll) ) has lead 
us to the cosmological model that is able to accommodate al- 
most all the av ailable high quality data- the so-called " concor- 
dance" model dBahcall et al.lll999USpergel et alJ2003h . Useful 
cosmological information has also been obtained from the Ly- 
a forest, weak lensing, galaxy cluster, and large-scale peculiar 
velocity studies. It is remarkable that this diversity of obser- 
vational data can be fully explained by a cosmological model 
that in its simplest form ha s only 5-6 free parameters (Liddle 
2004} iTegmark etaLll2004 . As the future data sets will be or- 
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ders of magnitude larger, leading to the extremely small statisti- 
cal errors, any further progress is possible only in case we fully 
understand various systematic uncertainties that could poten- 
tially bias our conclusions about the underlying cosmology. As 
such, one should try to use observables that are least sensitive 
to the theoretical uncertainties, contaminating foregrounds etc. 
Currently the "cleanest" constraints on cosmological models 
are provided by the measurements of the angular power spec- 
trum of the CMB. Since the underlying linear physics is well 
understood (see e.g. I Hull 1 995c iDodelso n 2003) we have a good 
knowledge of how the angular position and amplitude ratios 
of the acoustic peaks depend on various cosmological param- 
eters. However, the CMB data alone is able to provide accu- 
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rate measurements of only a few combinations of the cosmo- 
logical parameters. In order to break the degeneracies between 
the parameters one has to complement the CMB data with ad- 
ditional information from other independent sources e.g. the 
data from the type la supernovae, large-scale structure, or the 
Hubble parameter measurements. In fact, the well understood 
physical processes responsible for the prominent peak struc- 
ture in the CMB angular power spectrum are also predicted to 
leave imprints on the large-scale matter distribution. Recently 
the analysis of the spatial two-point correlation function of 
the Sloan Digit al Sky Survey (SPSS) Luminous Red Galaxy 
(LRG) sample (Eisenstein et al. 2005), and power spectra of 
the 2dF dCole et alJl2005h and SDSS LRG (JHUtsi 2006) red- 
shift samples, have lead to the detection of these acoustic fea- 
tures, providing the clearest support for the gravitational insta- 
bility picture, where the large-scale structure of the Universe is 
believed to arise through the gravitational amplification of the 
density fluctuations laid down in the very early Universe. 

In the current paper we work out the constraints on cos- 
mological par ameters using the SDSS LRG power spectrum as 
determined by Hiitsi (2006). In order to break the degeneracies 
between the parameters we complement our analysis with the 
data from other cosmological sources: the CMB data from the 
Wmap, and the measurement of the Hubble parameter by the 
HST Key Project 1 13 i We focus our attention on simple models 
with Gaussian adiabatic initial conditions. In the initial phase of 
the analysis we further assume spatial flatness, and also negligi- 
ble massive neutrino and gravitational wave contributions. This 
leads us to the models with 6 free parameters: total matter and 
baryonic matter density parameters: Q,„ and Q/,, the Hubble 
parameter h, the optical depth to the last-scattering surface t, 
the amplitude A s and spectral index n s of the scalar perturba- 
tion spectrum. |^j This minimal set is extended with the con- 
stant dark energy effective equation of state parameter w e s. We 
carry out our analysis in two parts. In the first part we use only 
the measurement of the ac oustic scale f rom the SDSS LRG 
power spectrum as given in iHiitsil (120061) . The analysis in the 
second part uses the full power spectrum measu r ement along 
with the covariance matrix as provided by Hiitsil (12006b . Here 
we add two extra parameters: bias parameter b and parameter 
Q that describes the deformation of the linear power spectrum 
to the nonlinear redshift-space spectrum. These extra parame- 
ters are treated as nuissance parameters and are marginalized 
over. Thus the larg est parameter space we should cope with is 
9-dimensional. [|£| Since the parameter space is relatively high 
dimensional it is natural to use Markov Chain Monte Carlo 
(MCMC) techniques. For this purpose we us e publicly avail- 
able cosmological MCMC engine Cosmomco ([Lewis & Bridle 
20021) . All the CMB spectra and matter transfer functions are 



calcu lated using the fast Boltzmann code Camh 17 ! (Lewis et al. 
20001) . 



13 http://www.ipac.caltech.edu/HOkp/ 

14 In fact, one might even consider a si mpler case with only 5 free 
parameters by fixing n s = 1 (Liddle 20041). 

15 Sinc e marginalization o ver the bias parameter can be done ana- 
lytically JBridle et al.ll2002l) the actual number of parameters can be 
reduced to 8. 

16 http://cosmologist.info/cosmomc/ 



The paper is organized as follows. In Sect. 2 we describe the 
observational data used for the parameter estimation. Sect. 3 
discusses and tests the accuracy of the transformations needed 
to convert the linear input spectrum to the observed redshift- 
space galaxy power spectrum. In Sect. 4 we present the main 
results of the cosmological parameter study and we conclude 
in Sect. 5. 

2. Data 

The SDSS LRG power spectrum as determined by Hiitsi ( 2006) 
is shown with filled circles and heavy solid errorbars in Fig.Q] 
There the upper data points correspond to the deconvolved ver- 
sion of the spectrum. Lij The thin solid lines represent the best- 
fitting model spectra, with the lower curve corresponding to 
the convolved case. As the survey window is relatively narrow 
the deconvolution can be done rather "cleanly". This decon- 
volved spectrum might be useful for the extra-fast parameter 
estimation employing analytic approximations for the matter 
transf er functions (Eisenstein & Hu 1998; iNovosyadlyj et al. 



1 19991) and fa st CMB angular powe r spectru m generators such 
as CMBfiH Sandvik et alJ l2004. D ASrMlKaplinghat et al 



2002b and CMBwarpFl djimenez et al.ll2004 . However, in this 
paper, as we use an accurate Boltzmann solver Camb to cal- 
culate CMB power spectra and matter transfer functions, the 
relative time taken by an extra convolution step is completely 
negligible. Thus in the following we use only the convolved 
spectrum. □ Accurate analytic fitting formulae for the sur- 
vey window functions can be found in Hiitsi (2006)0 The 
power spectrum covariance matrix in IHiitsil (120061) was mea- 
sured from 1000 mock catalogs generated with the second- 
order optimized Lagrangian perturbation calculation. The same 
paper also provides the measurement of the acoustic scale: 
(105.4 ± 2.3) hr l Mpc. This corresponds to the case when only 
sinusoidal modulation, as expected in the case of adiabatic ini- 
tial conditions, in the power spectrum is allowed. Relaxing this 
assumption by allowing an arbitrary phase shifts gave the re- 
sult, (103.0 ± 7.6) h Mpc, instead. In the followi ng parameter 
estimation process we use both of these values. In Hiitsi (2006) 
the measurement of the acoustic scale was achieved by first re- 
moving the "smooth" component of the spectrum and then fit- 
ting the parametrized family of functions to the oscillatory part 
via the modified version of the Levenberg-Marquardt method. 
The separation of the "smooth" and "oscillatory" components 
of the spectrum can be done rather accurately since the char- 
acteristic scales over which they change differ strongly. The 



17 http://camb.info/ 

1 8 The deconv olution was performed using an iterative algorithm due 
to Lucvl i 19741) with a specific implementation as given in lLin et alj 
(1996). 

19 http://www.hep.upenn.edu/ sandvik/CMBfit.html 

20 http://bubba.ucdavis.edu/DASh/ 

21 http://www.physics.upenn.edu/ raulj/CMBwarp/ 

22 Often also called a pseudospectrum. 

23 There the combination 'mode coupling kernels' is used in place 
of the more common 'window functions'. 
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Levenberg-Marquardt method which was used to determine 
the oscillation frequency approximates the likelihood surface 
near its maximum with a multidimensional Gaussian, and this 
way provides an approximate parameter covariance matrix. To 
avoid this "Gaussianity assumption" we have also performed 
a MCMC parameter estimation exercise, finding the best fit- 
ting acoustic scale along with its uncertainty in full agreement 
with the values quoted above. The question that might arise of 
course is how adequate is the parametric family that was used 
for fitting the oscillatory component? Even in the simplest case 
of the adiabatic initial fluctuations the damped sinusoidal mod- 
ulation is only an approximation. We investigate the possible 
biases introduced by assuming a fixed parametric form for the 
oscillatory part of the spectrum in more detail in Sect. 14.21 

As mentioned in the Introduction, in order to break several 
degeneracies between the cosmological parameters, we com- 
plement the SDSS LRG power spectrum data with the data 
from the Wmap CMB measurements. Specifica lly, we use the 



CMB temperature power spectrum as found in Hinsh aw et al 



(12003b a nd the tempe r ature- polarization cross-power as deter- 
mined bv lKogut et alj (12003b . The des cription of the likeli hood 
calculation using this data is given in IVerde et al.1 (|2003). We 
use the Fortran90 version of this likelihood code as provided 
by the Cosmomc package. 

While investigating the constraints arising from the mea- 
surement of the acoustic scale we do not run each time the 
full new MCMC calculation. Instead we importance sample 
the chains built for the Wmap data along with the constraint 
on the Hubble parameter as provided by the HST Key Project, 



H - 72 ± 8km/s/Mpc ( Freedman et alj 1200 lb . Using the 



Wmap data alone would result in too loose constraints on sev- 
eral parameters, and thus after importance sampling a large 
fraction of the chain elements would get negligible statistical 
weight, leaving us with too small effective number of samples. 

3. Power spectrum / acoustic scale transformation 

In this section we discuss the relation of the observed galaxy 
power spectrum to the underlying spectrum of the matter dis- 
tribution. We stress the need to take into account the so-called 
cosmological distortiono which almost always is being com- 
pletely neglected. q| This is fine for the very shallow surveys, 
but as we show later, for the samples like the SDSS LRGs, with 
an effective depth of Zes ~ 0.35, the cosmological distortion 
should certainly be taken into account. This is especially im- 
portant if power spectrum, instead of being well approximated 
by a simple power law, contains some characteristic features. 

There are other difficulties one has to face while trying to 
make cosmological inferences using the observed galaxy sam- 
ples. It is well known that galaxies need not faithfully fol- 



24 In order to convert the observed redshifts to the comoving dis- 
tances needed for the estimation of the power spectrum, one has to 
assume some background cosmological model - the fiducial cosmol- 
ogy. If the true underlying cosmology differs from the fiducial model 
we are left with a distortion of the power spectrum, which is often 
called the cosmological distortion. 

25 Acc ording to our kn owledge the only counter-example being the 
work bv lEisenstein et alj J2005 
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Fig. 1. Upper panel: Power spectra in somewhat unconven- 
tional form. Here the spectra have been multiplied by an ex- 
tra factor of k to increase the visibility of details. Filled circles 
with solid errorb ars rep r esent the SDSS LRG power spectrum 
as determined by iHiitsil (12006b . The upper data points provide 
the deconvolved version of the spectrum. The thin solid lines 
show the best-fitting model spectra. Lower panel: The same 
spectra as above now plotted in the usual form. 



low the underlyin g matte r distribution. This phenomenon is 






known as biasing (Kaiserl ll984l) . Whereas on the largest scales 
one might expect linear and scale -independent biasing (e.g. 
Coleslll993t iNaravanan et al.1 1200 0). on smaller scales this is 
definitely not the case. In general th e biasing can be scale - 
dependent, nonlinear, and stochastic (IDekel & Lahavl 1 1999b . 
The other complications involved are the redshift-space dis- 
tortions and the effects due to nonlinear evolution of the den- 
sity field. The redshift-space distortions, biasing, and nonlin- 
earities can be approximately treated in the framework of the 
Halo Model approach as described in Appendix lAl The imple- 
mentation of the Halo Model as presented there introduces four 
new parameters: M\ ov/ , the lower cutoff of the halo mass i.e. be- 
low that mass halos are assumed to be "dark"; a and Mo, the 
parameters describing the mean of the halo occupation distri- 
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bution i.e. the average number of galaxies per halo with mass 
M, which was assumed to have the form (N\M) - (jf) ; y, 
the parameter describing the amplitude of the virial motions 
inside the haloes with respect to the isothermal sphere model. 
This formulation of the Halo Mod el, along with the a ssumption 
of the best-fit Wmap cosmology (ISpergel et al.ll2003l) . is able to 
produce a very good fit to the observed SDSS LRG power spec- 
trum as demonstrated in Fig. Q] Moreover, all the parameters: 
Mi ow , a, Mo, y, are reasonably well determined. It turns out 
that to a good approximation these four extra parameters can 
be compressed down just to a single parameter Q, describing 
the deformation of the linearly evolved spectrum: 



Pgfiik) = b\\ + Qk?)Pm(k) . 



(1) 



Here b is the bias parameter and a good value for 77 turns out to 
be |. A similar type of parametric description for th e galaxy 



power spectrum was also used in ICole et al.l d2005l) . with a 



slight difference for the treatment of the largest scales. In that 
paper the authors suggest to take 77 = 2. However, we have 
found that 77 = | provides a better approximation for these 
4-parameter Halo Model spectra. |^| This is demonstrated in 
the upper panels of Figs. [2] and [3] There we have calculated 
a full range of Halo Model spectra (assuming the Wmap "con- 
cordance" cosmology) for different values of M\ ov , a, Mo, and 
y, drawn from the multidimensional Gaussian centered around 
the best-fitting values and with the parameter covariance ma- 
trix as found in Hiitsi (2006). Each of the calculated models is 
fitted with a simple parametric form as given in Eq. (Q]). The 
upper panels of Figs. [2]and[3]represent the density plots for the 
fractional accuracy of these simple fits i.e. for each wavenum- 
ber k they show the probability distribution functions for the 
achieved relative accuracy. With the heavy dashed lines we 
have also marked the 5% and 95% quantiles of the accuracy 
distributions. It is evident that r\ = | provides significantly 
better approximation to the spectral deformation than 77 — 2. 
The largest errors are seen to be located at the positions of the 
acoustic features, with a simple approximation in Eq. (Q} giving 
larger oscillation amplitudes. The Halo Model gives lower os- 
cillation amplitudes since relatively flat contribution from the 
1-halo term added to the 2-halo part starts to decrease the con- 
trast of the acoustic features, whereas the multiplicative trans- 
form in Eq. (Q]) preserves the contrast level of these wiggles. 
In the lower panels of Figs. [2] and [3] we have provided some 
examples of the pairs of spectra that correspond (staring from 
below) to the best matching case, and also the ones represent- 
ing the 68% and 90% quantiles of the distribution of the x 1 
values. The solid lines here correspond to the Halo Model cal- 
culations. For comparison also the SDSS LRG power spectrum 
along with the cubic spline fitted to the data points are shown. 
For clarity slight vertical shifts have been applied to the model 
spectra. As can be seen, the approximation in Eq. (Q]) is rather 
acceptable in the light of the accuracy of the SDSS LRG power 
spectrum measurement. This approximation is used in Sect. 14. 31 
where we fit the model spectra to the SDSS LRG data. 
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26 At least if the spectra have shapes close to the observed SDSS 
LRG spectrum. 



Fig. 2. Upper panel: A density plot showing the probability dis- 
tribution functions for the relative accuracy of the approxima- 
tion given in Eq. (Q]i with 77 = |. The set of Halo Model pa- 
rameters M] ow , a, Mo, and y, needed to calculate the "exact" 
spectra, were drawn from the multidimensional Gaussian dis- 
tribution centered at the best-fi t values and with a covariance 
matrix as found in iHiitsil (I2006I) . The heavy dashed lines mark 
the 5% and 95% quantiles of the relative accuracy distribu- 
tions. Lower panel: Filled circles with solid errorbars provide 
the SDSS LRG power spectrum. The data points are connected 
with a smooth cubic spline fit. The other set of lines represents 
some examples of the pairs of spectra that correspond, starting 
from below, to the best matching case, to the 68%, and to the 
90% quantiles of the distribution of the x 2 values. The solid 
lines show the Halo Model spectra while the dashed ones are 
the approximations from Eq. |[T). 



The cosmological distortion, mentioned in the beginning 
of this section, arises due to the simple fact that conversion of 
the observed redshifts to comoving distances requires the spec- 
ification of the cosmological model. If this cosmology differs 
from the true one, we are left with additional distortion of dis- 
tances along and perpendicular to the line of sight. In general, 
the spatial power spectrum measurements, in contrast to the 
angular spectra, are model dependent i.e. along with the mea- 
surements of the 3D power spectrum one always has to spec- 
ify the so-called fiducial model used to analyze the data. The 
fiducial model corresponding to the da ta shown in Fi g.[T]is the 
best-fit Wmap "concordance" model ( JSpergel et al.l 120031) . In 
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principle, for each of the fitted cosmological model one should 
redo the full power spectrum analysis to accommodate differ- 
ent distance-redshift relation. However, there is an easier way 
around: one can find an approximate analytical transformation 
that describes how the model spectrum should look like un- 
der the distance-redshift relation given by the fiducial model 
i.e. instead of transforming the data points we transform the 
fitted model spectra. Since the distance intervals along and per- 
pendicular to the line of sight transform differently, the initial 
isotropic theoretical spectrum P transforms to the 2D spectrum: 



P 2D (k h k ± ;z) 



1 



qi(z) • c{(z) 



\ c \\(z) 



k±_ 

cAz) 



;z 



(2) 



where the distortion parameters along and perpendicular to the 
line of sight are defined as: 



c ll(z) 
cjz) 



H na (z) 
H(z) 

djz) 

4 d (z) ' 



(3) 
(4) 



Here H(z) is the Hubble parameter and d ± (z) is the comoving 
angular diameter distance corresponding to the fitted theoreti- 
cal model. Superscript fid refers to the fiducial model. Here and 
in the following we use a tilde on top of P to denote theoretical 
spectrum "transformed to the reference frame of the fiducial 



cosmology". As we use the spectra that have the dimensions 
of volume an extra division by cy(z) • c]_(z) occurs due to the 
transformation of the volume elements: 



dV(z) = cn(z) • c\{z) ■ dV M (z) . 



By introducing the variables 

k J 
k 



*=V*ii + *i' ^ = t- 



»(z) 



c\\(z) 
cAz) 



(5) 



(6) 



we can express P 2D as follows: 



P 2D (k,v,z) 

1 



C||(z) ■ c 2 ± (z) 



c±(z) 



1 + 



1 



K 2 (z) - 1 



yu 2 ;z 



(7) 



Now the corresponding isotropized spectrum can be given as: 



P(k;z) 



1 



2q|(z)-cl(z) 



1 

I- 




dfi . (8) 



As the observations are done along the light-cone we have to 
perform relevant integrals along the redshift . The full treatment 
for th e light-cone effect can b e found in Yamamoto & Suto 
(11999l) : lYamamoto et al.l dl999l) . As we are investigating a two- 
point function, an accurate light-cone calculation would intro- 
duce two integrals ove r the redshifts (Ma tarrese et al.lll997at 
Yamamoto et al.l 1 1999b . However, it turns out that to a good 
approximation, excluding the very largest scales, the contribu- 
tions from different redshifts decouple and the double integral 
reduces to a simple one-dimensional integral over redshift. The 
final result for the P(k; z), averaged over the light-cone can be 
given as: 

/ ^dz ■ W 2 (k; z) n 2 (z) P(k; z) c\\(z) c\{z) 
P(k) = z -^- . (9) 

<,max 

J ^Idz • <W\k; z) h\z) c,|(z) c\{z) 



Here the result of lYamamoto et al.l (119991) has been generalized 
to include other weight factors in addition to the simple number 
density weighting. The most common weight functions 'Wiz) 
are the following: 



OV(/fc;z)oc 



i 
Hz) 
const 
l 



\+n(z)P(lc;z) 



for volume weighting 
for number weighting 
for the FKP weighting . 



(10) 



Here FKP stands for the weighting scheme due to 
Feldman et al. ( 1994). The power spectrum measurement of the 
SDSS LRGs in JHutsill2006l) used the FKP weighting function. 
In Eqs. (0 and ( TTOb h(z) represents the mean number density 
of galaxies as a funct i on of redshift. For the SDSS LRG sam- 
ple analyzed in Hiitsi (2006) the limiting redshifts z m ; n = 0.16 
and Zmax = 0.47. If instead of the integral over z in Eq. (0 we 
just take the integrand at the effective redshift (e.g. the median 
redshift) of the survey, and replace the distortion parameters c\\ 
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and c, with a single "isotropized" dilation of scales (see e.g. 
Eisenstein et al]|2005l) : 



Cisotr = y C||(Zeff) C^foff) , 

we can write instead of Eq. (0 simply 



1 / k 
P(k) = P[ 



(ID 



(12) 



Here the prefactor l/cj sotr can also be dropped, as it can be ab- 
sorbed into the bias parameter that is assumed to be a com- 
pletely free parameter throughout this paper. Although the true 
transformation for the power spectrum is different along and 
perpendicular to the line of sight, and also is dependent on red- 
shift, it turns out that a single dilation approximation taken at 
the median redshift of the survey can provide a very good ap- 
proximation, especially for relatively shallow surveys. F or the 
medi an redshift of the SDSS LRG sample as analyzed in Hiitsi 
(2006), z ~ 0.35, this approximation is very accurate as can 
be seen in Fig. [4] The upper panel of Fig. [4] shows a similar 
density plot as in Figs. [2] and [3] Here, in comparison to Figs. [2] 
and[3]where the background cosmology was fixed to the best-fit 
Wmap model and the Halo Model parameters were varied, we 
use the simple linear spectra while changing the cosmology. 
The set of cosmological models is drawn from the combined 
posterior corresponding to the Wmap plus HST Key Project 
data. As can be seen from this figure, for relatively shallow 
surveys the single "isotropized" dilation approximation is very 
precise: for ~ 90% of the models the approximation in Eq. 
( fT2l is more accurate than 0.5%. This is even more clear when 
looking at the lower panel of Fig. |4] where we have plotted the 
pairs of spectra corresponding to the best matching case, and 
also some examples representing 68% and 90% quantiles of 
the distribution of the x 2 values. As can be seen, even the pair 
of curves corresponding to the 90% quantile, are basically in- 
distinguishable. In Fig. [5] we have illustrated the case when the 
cosmological distortion is ignored. One can see that for ~ 90% 
of cases we make relative errors of ~ 6%, which is compa- 
rable to the amplitude of the acoustic oscillatory features. The 
lower panel of Fig.|5]presents pairs of spectra for 68% and 90% 
quantiles of the x 1 values. The inset shows the probability dis- 
tribution function for the "isotropized" dilation scale, as given 
in Eq. (fTTT i. compatible with the Wmap plus HST Key project 
constraints. Since the values of Ci sotl are quite often seen to dif- 
fer from Ci sotl - = 1 by 5 - 10%, it is clear that the cosmological 
distortion has to be taken into account if the power spectrum is 
measured as accurately as given by the SDSS LRG data points 
in the lower panel of Fig. [3] 

4. Results 

4. 1. Wmap + HST data 

As a starting point for several subsequent calculations we build 
a Markov cha in using the Wmap temperature-temperature 




0.13 oz OJS 

kpiMpa 1 ] 




Fig. 4. An analog of Fig [21 here provided in the context of 
the accuracy test for the cosmological distortion approximation 
given in Eq. d!21 l. The set of cosmological models was drawn 
from the combined posterior corresponding to the Wmap plus 
HST Key Project data. 



(Hinshaw et al.1 1 20031) and temperature-polarization 
(Kog ut et al.1 [2003) angular spectra in combination with 
the constraint on the Hubble parameter from the HST 



Key Project (Freedman et al.l 120011) . The results for the 2D 
marginalized distributions for all of the involved parameter 
pairs are shown in Fig. [6] Here the 68% and 95% credible 
regions are shown by solid lines. The original MCMC calcula- 
tion as performed by the Cosmomc software uses the variable 6 
- the angle subtended by the sound horizon at last scattering- 
in place of the more common Hubble parameter Hq. This 
leads to the better mixing of the resulti ng chain since 9 is onl y 
weakly correlated with other variables (lKosowskv et al.ll2.002h . 
The proposal distribution for all of the MCMC calculations 
carried out in this paper is taken to be a multivariate Gaussian. 
For the current Wmap + HST case we have used the CMB 
parameter covariance matrix as provided by the Cosmomc 
package. All of the seven default parameters here get implicit 
flat priors. The marginalized distributions in Fig.[6]are derived 
from a 100, 000-element Markov chain. As there is a very good 
proposal distribution available the chains typically equilibrate 
very fast and only a few hundred first elements need to be 
removed to eliminate the effects of the initial transients. We 
determine the length of this so-called burn-in p eriod using 
the Gibbsit software (iRaftery & Lewisl 1 19951) . The same 



http://www.stat.washington.edu/raftery/software.html 
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Fig. 5. As Fig. @] here instead showing the error one makes if 
cosmological distortion is completely neglected. In the lower 
panel we have shown only the examples corresponding to the 
68% and 90% quantiles. The inset shows the probability dis- 
tribution function for the "isotropized" dilation scale, as given 
in Eq. (fTTT i. compatible with the Wmap plus HST Key project 
constraints. 



program can also be used to estimate the length of the Markov 
chain required to achieve a desired accuracy for the parameter 
measurements. As a test one can run initially a short chain 
of a few thousand elements and analyze it with Gibbsit. It 
turned out that in the current case if we would like to achieve a 
1.25% accuracy at 95% confidence level for the measurement 
of the 2.5% and 97 .5% -quantiles of the most poorly sampled 
parameter, we would need a chain of ~ 25,000 elements. 
Thus according to this result our 100,000 element chain 
is certainly more than sufficient. Of course, all the various 
tools for diagnosing the convergence and for estimating the 
required chain length £j are just some more or less justified 
"recipes" that can lead to strongly incorrect results, especially 
in cases of poorly designed proposal distributions. Luckily, 
in cosmology as we have a very good knowledge about the 
possible parameter degeneracies, and also as the parameter 
spaces are relatively low dimensional, the construction of very 
good samplers is not too difficult. 



28 For a lot of online material related to these issues see 
http://www.statslab.cam.ac.uk/~mcmc/ . 
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Fig. 6. The 2D marginalized distributions for the Wmap + HST 
data. 

In the following subsections we use this Wmap + HST chain 
for the very fast determination of the parameter constraints re- 
sulting from the additional measurement of the SDSS LRG 
acoustic scale. The same chain was also used to produce Figs. 
E]and|5] 

4.2. Constraints from the measurement of the acoustic 
scale 

The low redshift acoustic scale as measured via the anal- 
ysis of the SDSS LRG power spectrum was found to be 
(105.4 ± 2.3) /i~' Mpc if adiabatic initial conditions were as- 
sumed (i.e. allowing only for the sinusoidal modulation in the 
spectrum), and (103.0 + 7.6) h~ l Mpc if this assumption was 
relaxe d by allowing additional oscillation phase shifts dHtitsi 
2006 ). These measureme nts refer to the Wmap best-fit cosmol- 
ogy dSpergel et al.ll2003l) which was used to analyze the SDSS 
LRG data. In Sect. [3] we have described accurate transfor- 
mations needed to accommodate other background cosmolo- 
gies. In the following we use SHI and SH2 to denote the 
sound horizon measurements of (105.4 + 2.3) /i -1 Mpc and 
(103.0 ± 7.6) h~ l Mpc, respectively. 

In this section we investigate the constraints on cosmo- 
logical parameters using the above given values for the sound 
horizon in combination with the Wmap data. We obtain initial 
bounds on parameters in a numerically efficient way by apply- 
ing the method of importance sampling on the earlier calcu- 
lated Wmap + HST chain. However, to be confident in the re- 
sults obtained we always carry out a full MCMC calculation 
from scratch for each of the considered cases. As a final re- 
sults we only quote the constraints on cosmological parameters 
obtained from the direct MCMC calculations. Importance sam- 
pling is only used as an independent check of the validity of the 
results. In general both methods reach to the parameter bounds 
that are in a good agreement. 



G. Hiitsi: SDSS LRG cosmological parameters 





















- 


-~T 


''■Vp"^) 


tij/*?*?'**** 




"**: 


-4A- 


;^~i ? 


-rr- 


K- 


:/..".v:' v ,' 


- 


-0.77 h 


Mpc 














- 





400 600 

sample # 



0.004 
0.003 
0.002 


I I I I I 




II 

' \ ' \ ' — 




1 


■■■ 


0.001 


1 


/ 




^ 


- 


■/\ 





' - 


N / "\ 


0.001 
0.002 
0.003 


/ 


/ 


V 


0.004 


N^^Hl' 




1 1 1 1 1 



1700 
600 
500 
400 
- H 300 
200 
100 




85 90 95 100 105 110 115 120 125 130 
s«ed[h" 1 Mpc] 



Fig. 7. Upper panel: Comparison of the sound horizon as de- 
termined from ~ 1000 mo del spectra, v ia the same fitting tech- 
niques that were used in iHiitsil (120061) to measure the SDSS 
LRG sound horizon, with the analytical approximation given 
in Eqs. ED, dOt , dB~9l , dBTTOb . ( fB~TTI The model spectra 
were drawn from the posterior distribution corresponding to 
the Wmap + HST data. Lower panel: The density plot of the 
residuals after removing the average bias of 0.77 h Mpc. The 
solid dashed lines mark the 68% credible region. 



It is fine to use importance sampling if new constraints are 
not too constraining and are consistent with the earlier gener- 
ated chain. Having a measurement of the acoustic scale TPl 
with an error AT, importance sampling simply amounts to mul- 
tiplying each original sample weight by 



fi = exp 



v^model/ $) 



2 AT 



(13) 



where Tmodeij denotes the theoretical sound horizon correspond- 
ing to the i-th Markov chain element. The physical size of the 
sound horizon s at the end of the drag-epoch is is determined 
by the parameter combinations D.„,h 2 and Q.bh 2 i.e. physical 
densities of the CDM and baryonic c omponents. Accurate fit - 
ting formulae fo r s can be found in iHu & Sugiyamal (11996); 
Eisenstein & Hul Ql_998). We have provided these in Appendix 



|B1 where also the transformation into different cosmological 



29 We use tilde to denote the quantities that are "tied to the" fiducial 
cosmological model used to analyze the data. 



frame is described. This transformation induces an extra de- 
pendence of the sound horizon T as measured from the matter 
power spectrum, on h and w e ff. For more details see Appendix 
IB1 The dependence of Tat redshift z ~ 0.35 on various parame- 
ters for spatially flat models around the best fitting Wmap model 
point can be conveniently expressed as the following principal 
component: 



0.14 



-0.28 



0.022 



-o.io 



h 
071 



0.94 , 



014 AT 

= 1 ± — . (14) 



As probably expected, for those relatively small redshifts by 
far the strongest dependence is on the Hubble parameter h. 
To avoid any biases due to the approximate nature of the Eqs. 
( lB.9t , dB~T0l i.( lB. 1 ll ). and also due to the specific method used 
to measure the oscillation frequency in the SDSS LRG power 
spectrum, we carry out the following Monte Carlo study. We 
draw ~ 1000 samples from the Wmap + HST chain by thin- 
ning it by a factor of ~ 10. For each of the parameter combina- 
tions we calculate theoretical matter spectra using Camb. The 
oscillatory components of the spectra are extracted by dividing 
th em with a "smoo thed" a ppr oximate model spectra as given 
i lEisenstein & Huldl998l) . 



, £3 The resulting "flattened" spectra 
are fitted with damped sinusoidal waves |^j and the sound hori- 



zon Tfltted is determined via the Levenberg-Marquardt fitting 
technique. All the spectra are calculated at exactly the same 
wavenumbers as the data points given in Fig [TJ The power 
spectrum co variance matrix is taken from the Appendix G of 
Hiitsi (2006). For each model the sound horizon Tnodei is calcu- 
lated using Eqs. dB~4l (|B~B1 ( |B~9l , dBTOl i, dBTO . The com- 
parison of Tfitted versus T mo d e i is provided in Fig. [7] In the upper 
panel we have plotted Tfj tt ed - Tnodei- As can be seen there is 
a slight tendency for the fitted values T ntt ed to be larger than 
Tnodei- After removing the constant bias of 0.77 h Mpc the 
remaining fluctuations are < 0.5%, which is demonstrated in 
the lower panel of Fig. [7] This plot is an analog to the earlier 
density plots shown in Figs. [2 [3j IH |5] Here the dashed lines 
show the region covering 68% of all the cases. These were the 
results that apply to the case when the phase of the sinusoidal 
waves was not allowed to vary. If the phase is additionally al- 
lowed to change the corresponding effective bias turns out to 
be in the opposite direction with a value 1.8 h' 1 Mpc, instead. 
Thus the bias corrected values for the sound horizon used in 
our analysis are the following: (104.6 ± 2.3) /r 1 Mpc (SHI) 
and (104.8 ± 7.6) /T 1 Mpc (SH2). One might even go a step 
further and instead of removing only a constant offset, remove 
also the next order i.e. the linear component. This more accu- 
rate treatment has probably rather negligible effect on the fi- 
nal results, since around the measured sound horizon values of 
~ 105 /z" 1 Mpc the accuracy after removing the constant offset 
is already ~ 0.2-0.3%, which is an order of magnitude smaller 
than the measurement errors of 2.3 - 7.6 hr x Mpc. 



30 The separation of the oscillatory component and the underlying 
smooth CDM continuum can be done very cleanly due to significantly 
different characteristic scales over which they change. The small resid- 
ual deformations of the oscillatory part have negligible impact on the 
inferred oscillation period. 

31 For a precise parametric form see Hiitsi (2006). 
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Fig. 8. The 2D marginalized distributions for the Wmap data Fig. 9. The same as Fig. [8] only for the sound horizon 
along with the constraint on the low redshift sound horizon, ment (104.8 ± 7.6) h' 1 Mpc (SH2). 
J - (104.6 + 2.3) h' 1 Mpc (SHI), obtained via the importance 
sampling of the Wmap + HST results shown in Fig. [6] 



measure- 



Using this correction for the bias and the method to calcu- 
late the theoretical size of the sound horizon at the end of the 
drag-epoch, as presented in Appendix iBl we can immediately 
perform the relevant reweighting of the Wmap + HST chain 
(see Eq. (TT3I )). It turns out that relatively large fraction of the 
Wmap + HST chain elements "survive" this reweighting proce- 
dure, justifying the use of the importance sampling method. In 
particular, for the SH2 case we are left with ~ 36, 000, and for 
the SHI ~ 12, 000 samples. The results of this calculation in the 
form of the 2D marginalized distributions is presented in Figs. 
[8]and|9] Here Fig. ISI91 corresponds to the SH1/SH2 case. In 
comparison to the analogous Fig.|6]the most dramatic changes 
are for Hq and w e g, whereas the rest of the parameters stay es- 
sentially the same. The HST constraint for the initial Wmap 
chain was just implemented in order not to loose too many 
samples in current importance sampling calculations. As can 
be seen, the new constraints on Hq are significantly stronger 
than the one provided by the HST. In Fig. [8] due to somewhat 
lower number of samples (~ 12, 000) the contours start to be- 
come more noisy. Earlier we estimated that the measurement 
of the 2.5% and 97.5% quantiles with an accuracy of 1.25% at 
95% CL requires ~ 25, 000 samples. So is this 12, 000 enough 
for the parameter estimation purposes? To test that we have 
also performed a full MCMC calculation from scratch (with 
100, 000 samples) using Wmap data along with a sound horizon 
measurement SHI. The results of this calculation are shown in 
Fig. [10] The contours in Fig. [8] although being noisier are very 
similar to the ones in Fig.QJJ In fact the corresponding ID dis- 
tributions are practically indistinguishable. This shows that the 
initial use of importance sampling was indeed justified. 
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Fig. 10. The exact analog of Fig. g] now for the full MCMC 
calculation. 



4.3. Constraints from the full power spectrum 

Using the Wmap data and the SDSS LRG power spectrum as 
shown in Fig. [TJ along with the power spectrum transformation 
and an additional new parameter Q, as described in Sect. [3] we 
build a 100,000 element Markov chain in the 8-dimensional 
parameter space. The resulting 2D parameter distribution func- 
tions are shown in Fig.Qj] Here we see that in several cases dis- 
tributions start to become doubly-peaked. Also the constraints 
on Hq and w e g are weaker than the ones obtained in the pre- 
vious subsection. On the other hand, now a rather strong con- 
straint has been obtained for Q. c d m h 2 . Even stronger constraint 
(not shown in the figure) is obtained for Q,„/z- the shape pa- 
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Fig. 11. The 2D marginalized distributions from the Wmap + 
SDSS LRG full power spectrum MCMC calculation. 



rameter F. This just illustrates the the well-known fact that 
the shape of the matter power spectrum is most sensitive to T. 
The new parameter Q, describing the deformation of the linear 
spectrum to the evolved redshift-space galaxy power spectrum, 
is seen to be significantly degenerate with several parameters 
e.g. Qi/1 2 , n„ r, A s . On the other hand it does not interfere too 
strongly with Ho- 

It might seem strange that using the full data we obtain 
weaker constraints on Ho and w e ff. But after all, we should 
not be too surprised, since our understanding of how the lin- 
ear spectrum is deformed to the evolved redshift-space power 
spectrum is rather limited. Here we were introducing an addi- 
tional parameter Q, which starts to interfere with the rest of the 
parameter estimation. Also one should remind that maximum 
likelihood is the global fitting technique i.e. it is not very sensi- 
tive to specific features in the data. On the other hand, model- 
ing of the oscillatory component of the spectrum does not call 
for any extra parameters. Also the underlying physics is much 
better understood. In fact, the observable low redshift acoustic 
scale is determined by four parameters only: Q. m h 2 , Q/,/z 2 , w e ff 
and h. The optimal data analysis of course should incorporate 
both components: (i) general shape of the spectrum i.e. low fre- 
quency components, and (ii) oscillatory part, with appropriate 
weightings. It is clear that in the current "full spectrum" maxi- 
mum likelihood analysis the acoustic features are weighted too 
weakly. 

4.4. One dimensional distributions 

To compare the measurements of the parameters in a more clear 
fashion we provide in Fig.[T2"lseveral ID marginalized distribu- 
tions. The 68% and 95% credible regions along with the me- 
dians of these distributions are provided in Table Q] Here the 
parameters D./,h 2 , Q. c dmh 2 , 6, T, w e s, n s , A s and Q (the last in 
case of the full spectrum analysis only) are primary parameters 
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Fig. 12. The ID posterior distributions for several cosmological 
parameters. Solid, dash-dotted, dashed, and dotted lines corre- 
spond to the Wmap + HST, Wmap + SDSS LRG SHI, Wmap + 
SDSS LRG SH2, and Wmap + SDSS LRG full power spectrum 
cases, respectively. The compact summary of these results can 
be found in Table Q] 

as used in the MCMC calculations. All the rest: Qa, to, Q m , 
Zreion, Hq, qo, jo are derived from these. Here to is the age of the 
Universe, qp the deceleration par ameter and jo the so-called 
jerk (see e.g. lBlandford et al.ll2005l) at z - 0. The deceleration 
parameter qo and jerk jo are introduced as usual via the Taylor 
expansion of the scale factor: 

-qoHl(t - to) 2 + 



a(t) = ao 



1 + H (t - t ) 



2 



fo) 3 + 



(15) 



From Fig.Q~2]we can see that many parameters stay essentially 
the same as determined by Wmap + HST data. On the other 
hand, a new precise measurement of Ho, thanks to the mea- 
surement of the low redshift sound horizon along with strong 
constraints on Q/,/i 2 and Q«/,„/z 2 from the CMB data, helps to 
determine Q,„ (as well as Q/, and Cicdm separately) rather pre- 
cisely. The same applies to the case of the full spectrum anal- 
ysis, which provides us with a good estimate for the shape pa- 
rameter T = Q. m h. In both cases also the constraint on w e g is 
significantly improved. New improved limits on Q,„ and w e g 
immediately transform to better constraints on qo and jo (see 
Appendix 0. For the "vanilla" ACDM model with w — -1 
the jerk parameter jo = 1. We can see that at the moment 
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jerk is still rather poorly constrained. Using only the observa- 
tional data whose nature is very well understood, namely the 
CMB power spectra along with the low redshift sound horizon 
measurement, we get very strong support for the accelerating 
Universe (i.e. qo < 0). The values qo > are ruled out by 1.4cr, 
2.9o- and 5.5cr in case of the Wmap + HST, Wmap + SDSS 
LRG SH2 and Wmap + SDSS LRG SHI, respectively. Of 
course, one has to remind that until now the analysis assumed 
flat spatial sections. 

4.5. Most interesting constraints 

We have shown that by adding the SDSS LRG clustering data 
to the Wmap results we can get significantly tighter constraints 
on //o, Q. m and w a s (or qo and jo) than from the Wmap + HST 
analysis alone. The comparison of the obtained limits on pa- 
rameters Ho, Q. m and w e ff is provided in Fig. [13] The largest 
error contours in both upper and lower group of panels corre- 
spond to Wmap + HST, while the tightest to the Wmap + SDSS 
LRG SHI case. In the upper group of panels we have addition- 
ally given the constraints for the Wmap + SDSS LRG SH2 case, 
whereas the lower group provides extra limits from the full 
spectrum + Wmap analysis. In addition, in each of the panels 
we have given the degeneracy lines corresponding to the princi- 
pal component given in Eq. (fT4l) .l 33 lln Q m - Hq plane we have 
additionally plotted the lines corresponding to Q„,/z 2 = const 
and T = Q. m h - const. These are the combinations well deter- 
mined by the CMB data and by the general shape of the matter 
power spectrum, respectively. As is evident from Fig. [13J the 
principal direction of the low redshift sound horizon constraint 
is always almost perpendicular to the corresponding Wmap + 
HST error contours, demonstrating the high level of comple- 
mentarity of this new measurement. For the spatially flat mod- 
els with constant dark energy equation of state parameter there 
exists a unique relation between parameter pairs (Q,„, w e s) and 
(qo, jo) (see AppendixO. Fig.[l4]presents similar plots to Fig. 
fT3l now only for the parameter triad (Ho, qo, jo) instead. The 
parameters shown in Figs.Q~3]and[T4]are the ones that determine 
the low redshift expansion law. Introducing the look-back time 
fib - to- t, where to is the age of the Universe at z - 0, one can 
find for the redshift: 

z - H t lh + (l + y ) #o«?b + (l + *> + f )#o*?b + ■ • ■ (16) 

The precise calculation for the look-back time as a function of 
redshift is shown in Fig. [15] Here the upper panel shows the 2<x 
regions corresponding to the Wmap + HST and Wmap + SDSS 
LRG SHI, respectively. The inset in the upper panel displays 
these regions after dividing by the look-back time correspond- 
ing to the best-fit Wmap cosmology. Here in addition to the 2<x 
contours also lcr regions are given. It is evident that the low 
redshift sound horizon measurement has helped to determine 
the recent expansion history of the Universe with much greater 
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32 We can perform this analysis of the far away tails of the distri- 
butions since the Wmap + HST chain we start with contains enough 
samples with q > (see Fig.ll2ll. 

33 The analog of Eq. d 14b valid for the non-flat cases is given in 
Appendix iBl 



60 70 80 90 0.1 0.2 0.3 0.4 0.5 0.6 



Fig. 13. The comparison of constraints on Ho, Q. m and w e fi. In 
all panels the largest error contours correspond to the Wmap + 
HST, while the tightest to the Wmap + SDSS LRG SHI case. 
The upper group of panels shows additionally the constraints 
for the Wmap + SDSS LRG SH2 case, whereas the lower group 
provides extra limits from the full spectrum + Wmap analysis. 
The dashed lines in all the panels show the principal component 
from Eq. ( [141 ). The additional lines in Q. m - Hq plane provide 
the directions Q m /z 2 = const and F = Q. m h - const. 



accuracy than available from the Wmap + HST data alone. Of 
course, this is largely due to the much tighter constraint ob- 
tained for the Hubble parameter. The lower panel in Fig. [15] 
shows a similar plot than the inset in the upper panel. Here 
we have given only the lcr regions as a function of redshift 
for (starting from the bottommost) the Wmap + HST, Wmap + 
SDSS LRG full spectrum, Wmap + SDSS LRG SH2 and Wmap 
+ SDSS LRG SHI cases. 
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Table 1. Various quantiles of the ID distributions shown in Fig. [12] The first group of parameters are the primary ones used in 
the MCMC calculations, the second group represents various derived quantities, and the last shows the parameters held fixed due 
to our prior assumptions. The last row of the table also gives the total number of free parameters, excluding the bias parameter 
that was marginalized out analytically, for all of the investigated cases. Also shown are the ^-values for the best-fitting model 
and the effective number of degrees of freedom involved. 



Parameter 


WMAP + HST 


WMAP + SDSS LRG 
sound horizon; adiab. (SHI) 


WMAP + SDSS LRG 
sound horizon (SH2) 


WMAP + SDSS LRG 
full spectrum 


Median 


68% 


95% 


Median 


68% 


95% 


Median 


68% 


95% 


Median 


68% 


95% 


n b h 2 

e 

T 

n s 

log(10 10 A,) 
1 


0.0238 

0.120 

1.0489 

0.153 

-0.97 

0.991 

3.24 


+0.0026 

-0.0014 

+0.017 

-0.018 

+0.0085 

-0.0073 

+0.044 

-0.153 

+0.30 

-0.36 

+0.084 

-0.037 

+0.28 

-0.16 


+0.0057 

-0.0025 

+0.031 

-0.033 

+0.0179 

-0.0138 

+0.231 

-0.153 

+0.54 

-0.73 

+0.183 

-0.061 

+0.53 

-0.28 


0.0239 

0.119 

1.0496 

0.160 

-0.96 

0.993 

3.26 


+0.0028 

-0.0014 

+0.018 

-0.018 

+0.0088 

-0.0081 

+0.047 

-0.160 

+0.18 

-0.22 

+0.093 

-0.037 

+0.31 

-0.18 


+0.0062 

-0.0027 

+0.033 

-0.035 

+0.0182 

-0.0149 

+0.238 

-0.160 

+0.32 

-0.45 

+0.195 

-0.063 

+0.53 

-0.32 


0.0238 

0.118 

1.0488 

0.161 

-0.95 

0.992 

3.24 


+0.0024 

-0.0014 

+0.018 

-0.017 

+0.0083 

-0.0070 

+0.043 

-0.161 

+0.20 

-0.23 

+0.083 

-0.036 

+0.28 

-0.17 


+0.061 

-0.025 

+0.036 

-0.033 

+0.0178 

-0.0131 

+0.239 

-0.161 

+0.34 

-0.53 

+0.190 

-0.061 

+0.54 

-0.29 


0.0237 

0.118 

1.0486 

0.152 

-0.86 

0.993 

3.22 

1.71 


+0.0025 

-0.0013 

+0.011 

-0.010 

+0.0075 

-0.0065 

+0.046 

-0.152 

+0.21 

-0.40 

+0.092 

-0.036 

+0.30 

-0.16 

+0.63 

-0.60 


+0.0043 

-0.0023 

+0.020 

-0.018 

+0.0146 

-0.0119 

+0.195 

-0.152 

+0.34 

-0.92 

+0.141 

-0.060 

+0.45 

-0.24 

+1.25 

-1.11 


tfo[km/s/Mpc] 
k 

1 L m 

to 

£reion 

r = Q, m h 


71.4 

-0.54 

0.91 

0.717 

0.283 

13.42 

16.5 

0.202 


+8.0 

-8.2 

+0.38 

-0.45 

+1.55 

-0.54 

+0.063 

-0.083 

+0.083 

-0.063 

+0.31 

-0.36 

+7.3 

-5.5 

+0.034 

-0.031 


+ 14.9 

-15.5 

+0.68 

-0.89 

+4.11 

-0.69 

+0.108 

-0.189 

+0.189 

-0.108 

+0.64 

-0.76 

+ 10.9 

-10.3 

+0.073 

-0.058 


70.8 

-0.53 

0.87 

0.715 

0.285 

13.38 

16.9 

0.202 


+2.1 

-2.0 

+0.16 

-0.18 

+0.76 

-0.45 

+0.031 

-0.032 

+0.032 

-0.030 

+0.34 

-0.40 

+7.8 

-5.7 

+0.024 

-0.023 


+4.4 

-4.0 

+0.29 

-0.37 

+1.82 

-0.67 

+0.062 

-0.059 

+0.059 

-0.062 

+0.65 

-0.84 

+ 10.6 

-10.7 

+0.044 

-0.044 


70.5 

-0.51 

0.85 

0.712 

0.288 

13.43 

16.8 

0.202 


+3.8 

-3.7 

+0.20 

-0.24 

+0.83 

-0.47 

+0.040 

-0.042 

+0.042 

-0.040 

+0.32 

-0.39 

+7.2 

-5.9 

+0.026 

-0.024 


+7.8 

-7.3 

+0.35 

-0.51 

+2.31 

-0.65 

+0.076 

-0.090 

+0.089 

-0.077 

+0.60 

-0.87 

+11.0 

-10.5 

+0.050 

-0.046 


67.6 

-0.38 

0.63 

0.693 

0.307 

13.48 

16.4 

0.207 


+7.7 

-4.3 

+0.24 

-0.50 

+1.43 

-0.33 

+0.045 

-0.034 

+0.034 

-0.045 

+0.27 

-0.26 

+7.7 

-5.9 

+0.011 

-0.0P 


+ 17.9 

-7.2 

+0.38 

-1.22 

+5.31 

-0.39 

+0.096 

-0.064 

+0.064 

-0.096 

+0.53 

-0.52 

+ 10.0 

-10.6 

+0.023 

-0.0^7 


n k 

Em„[eV] 
At 

dln«. T 
dint 










































































^ 2 /dof 


1429/1404 


1429/1404 


1429/1404 


1439/1417 




7 free parameters 


7 free parameters 


7 free parameters 


8 free parameters 



5. Extended analysis 

In this section we investigate the effect of several previously 
made assumptions and carry out an extended analysis by relax- 
ing some of these. 

We start out with a small comment about tensor modes. 
There exist several inflationary models which predict a non- 
negligible tensor mode contribution to the CMB angular 
power spectrum (for a classification of several models see e.g. 
Dodelson et ajj 1 19971) . However, the inclusion of the tensor 
modes into our analysis would not help to constrain them bet- 
ter than the results obtainable from the CMB data alone Q In 
order to make a clean separation of the tensor and scalar con- 
tributions to the CMB angular fluctuations one would greatly 
benefit from the independent knowledge of the scalar fluctu- 
ation level as obtainable from the large-scale structure stud- 
ies. Nevertheless, as in our analysis the biasing parameter is 
treated as a completely free quantity that is marginalized out, 
we do not have any sensitivity to the absolute level of the scalar 
fluctuation component. It is worth pointing out that in princi- 
ple a good handle on bias parameter can be obtained by the 
study of the higher order clustering meas ures e.g. bispectrum 



study or the higher order clustering measu 
dMatarrese et alll997bHvferde et alll998l) . 



34 Fo r the results using the CMB data alone see e.g. ISpergel et al.l 
J2003h . 



Since several neutrino oscillation experiments unambigu- 
ously indicate t hat neutrinos have a non-zero mass (for a recent 
review see e.g. Lesgourgues & Pastor 2006) it would be inter- 
esting to investigate models with massive neutrino component. 
For simplicity we concentrate on models with three generations 
of neutrinos with degenerate masses. Also, as the generic pre- 
diction of almost all the inflationary models is the nonmeasur- 
ably small spatial curvature i.e. Q.^ a* it is of great interest 
to test whether this is compatible with the observational data. 
For these reasons we have extended our initial analysis to allow 
for the massive neutrino component and also have investigated 
the effect of relaxing the assumption about spatial flatness. The 
results of this study for some of the model classes are briefly 
presented in Table [2] Here we have not carried out an analysis 
using the full measured LRG power spectrum, instead only the 
measurement of the acoustic scale was added as an additional 
information to the CMB data. This can be justified for two rea- 
sons: (i) To model the observed power spectrum correctly one 
needs to incorporate the corrections due to nonlinearities and 
redshift-space distortions, which in our case was done sim- 
ply by introducing one additional parameter Q (see Eq. (Q3). 
However, it is clear that the effect of this parameter can resem- 
ble very closely the power-damping effect of massive neutrinos 
i.e. there will be very strong degeneracies that result in poor 
constraints on neutrino mass. Instead of the dynamic damp- 
ing effect of massive neutrinos we can probe their effect on 
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Table 2. Constraints on selected parameters from an extended analysis. Also shown are the ^-values for the best-fitting model 
and the effective number of degrees of freedom involved. 





7 free parameters: H/,/r, Q. cc i m h z , 8, r, n s , A s , H,t 




fixed parameters: iVeff = — 1, £m r - 0, At - 0, ";'}„'/ = 


Parameter 


WMAP + HST 


WMAP + SDSS LRG 
sound horizon; adiab. (SHI) 


WMAP + SDSS LRG 
sound horizon (SH2) 


Median 


68% 


95% 


Median 


68% 


95% 


Median 


68% 


95% 


ff [km/s/Mpc] 

n,„ 
n t 

10 

r = n„,/i 


67.8 
0.261 
0.767 
-0.029 

-0.64 
0.179 


-7.9 

+0.107 

-0.068 

+0.074 

-0.110 

+0.023 

-0.025 

+0.16 

-0.10 

+0.051 

-0.038 


-15.1 

+0.246 

-0.114 

+0.119 

-0.226 

+0.041 

-0.060 

+0.36 

-0.17 

+0.098 

-0 065 


66.8 
0.261 

0.777 
-0.037 
-0.647 
0.174 


-28 

+0.027 
-0.028 
+0.049 
-0.056 
+0.025 
-0.024 
+0.069 

+8:8§1 

-0 076 


-6.9 

+0.048 

-0.053 

+0.090 

-0.092 

+0.047 

-0.042 

+0.114 

;8:oM 

-0.048 


65.9 

0.266 

0.768 

-0.036 

-0.635 

0.177 


-4.8 

+0.049 

-0.041 

+0.061 

-0.065 

+0.026 

-0.024 

+0.088 

-0.082 

+0.030 

-0.031 


-9.9 

+0.101 

-0.074 

+0.102 

-0.120 

+0.047 

-0.053 

+0.170 

-0.137 

+0.057 

-0 05? 


^/dof 


1429/1404 


1429/1404 


1429/1404 





8 free parameters: €lhh A , Q^^h 1 , 0, r, n s , A s , vv c |t, Q^- 




fixed parameters: Sm,, - 0, Aj - 0, a ^n$ = ^ 


Parameter 


WMAP + HST 


WMAP + SDSS LRG 
sound horizon; adiab. (SHI) 


WMAP + SDSS LRG 
sound horizon (SH2) 


Median 


68% 


95% 


Median 


68% 


95% 


Median 


68% 


95% 


fl [km/s/Mpc] 
"a 

10 

r = n„,/i 


63.3 

0.32 

0.71 

-0.94 

-0.023 

-0.51 

0.206 


M 

-0.46 
+0.87 
-0.61 
+0.025 

Ml 2 

-0.70 
+0.160 

-0050 


+18.4 

~M 
+8:!f 
+8f 

-0.99 
+0.088 
-0.071 
+1.01 

+olh 

-0 087 


65.6 
0.292 

0.744 

-1.20 

-0.037 

-0.83 

0.190 


+4.0 

-4-4 

+0.043 

-0.05 1 

+0.35 

-0.50 

+0.024 

-0.028 

+0.34 

-0.47 

+0.030 

-0 036 


+8.1 

+0:073 

+8:??7 

-0.079 

+0.51 

-0.76 

+0.048 

-0.065 

+0.52 

-0.80 

+0.049 

-0 065 


64.8 

0.291 

0.746 

-1.13 

-0.037 

-0.75 

0.190 


+5.4 

+o:06o 

+8:8S 

-0.061 

+0.31 

-0.54 

+0.024 

-0.035 

+0.32 

+8:oi2 

-0 014 


+11.4 
-11.4 
+0.118 

+8:??7 

-0.115 

+0.47 

-0.81 

+0.052 

-0.070 

+0.52 

-0.83 

+0.062 

-0 06? 


r/dof 


1428/1403 


1427/1403 


1427/1403 





7 free parameters: l\/r, Q. Ci i, u h 1 , 9, r, n s , A s , Sm v 




fixed parameters: w^\ = — 1, £1% = 0,At = 0, tl ""/ = 


Parameter 


WMAP + HST 


WMAP + SDSS LRG 
sound horizon; adiab. (SHI) 


WMAP + SDSS LRG 
sound horizon (SH2) 


Median 


68% 


95% 


Median 


68% 


95% 


Median 


68% 


95% 


ff [km/s/Mpc] 
Em„[eV] 

10 

r = n„,/i 


69.2 
0.300 

-0.54 

0.208 


-6.1 

+0.086 

-0.066 

< 0.80 

+0.13 
-0.10 
+0.037 
-0 034 


-11.2 
+0.174 
-0.118 
< 1.71 
+0.27 
-0.18 
+0.077 
-0 060 


69.2 
0.294 

-0.552 
0.204 


-2.5 

+0.014 

-0.013 

<0.72 

+0.024 
-0.021 
+0.009 
-0 009 


-5.4 

+0.028 

-0.025 

< 1.63 

+0.050 
-0.041 
+0.018 
-0 017 


68.5 
0.302 

-0.539 

0.207 


-3.9 

+0.044 

-0.039 

<0.77 

+0.069 

+8:8S1 

-0 019 


-7.5 

+0.094 

-0.072 

< 1.67 

+0.150 
-0.112 
+0.042 
-0 035 


*-/dof 


1429/1404 


1429/1404 


1429/1404 





8 free parameters: Cl/ ? h~, Q-dmh 1 , 0, t, n s , A s , w C |T, l,m v 




fixed parameters: Q.^ = 0, Aj - 0, ".,'""/ - 


Parameter 


WMAP + HST 


WMAP + SDSS LRG 
sound horizon; adiab. (SHI) 


WMAP + SDSS LRG 
sound horizon (SH2) 


Median 


68% 


95% 


Median 


68% 


95% 


Median 


68% 


95% 


H [km/s/Mpc] 

n,„ 

M-'elT 

Im,,[eV] 

10 

r = n,„h 


70.0 
0.298 

-1.14 

-0.68 
0.209 


-7.3 

+0.093 

-0.073 

+0.38 

-0.42 

<0.95 

+0.42 
-0.47 
+0.040 
-0 035 


-14.8 
+0.226 
-0.122 
+0.66 
-0.74 

< 1.80 

+0.80 
-0.92 
+0.090 
-0 061 


69.2 
0.306 

-1.08 

-0.62 
0.211 


-2.8 

+0.044 

-0.051 

+0.31 

-0.37 

< 0.95 

+0.26 
-0.28 
+0.027 
-0 01! 


-5.3 

+0.076 

-0.085 

+0.45 

-0.74 

< 1.80 

+0.39 
-0.56 
+0.046 
-0 05? 


68.6 
0.313 

-1.12 

-0.64 
0.214 


-3.6 

+0.049 

-0.052 

+0.31 

-0.38 

<0.95 

+0.28 
-0.32 
+0.028 
-0 037 


-7.2 

+0.100 

-0.095 

+0.48 

-0.72 

< 1.80 

+0.46 
-0.60 
+0.052 
-0 051 


;r/dof 


1429/1403 


1429/1403 


1429/1403 



the kinematics of the background expansion by exploiting the 
measurement of the low redshift acoustic scale, (ii) The non- 
zero spatial curvature has twofold effect on the matter power 
spectrum. First, it influences the growth rate of the fluctuations 
and thus changes the amplitude of the spectrum. Since in our 
case the amplitude is treated as a free parameter we are unable 
to use this effect. Second, the spatial scales get transformed, 
resulting in the horizontal stretch/compression of the power 
spectrum. It is clear that the smooth "continuum" of the power 
spectrum with an unknown amplitude does not contain much 
information about the possible horizontal stretch/compression 
l 35 l as this transformation can be easily mimicked by the cor- 



35 That is especially true in the case we do not have information 
about the scale of the turnover in the spectrum (see Fig.[T]lower panel). 



responding change of the normalization. The degeneracy be- 
tween these two transformations can be broken if the spectrum 
contains some sharper features e.g. acoustic oscillations. 



According to the results presented in Fig. IB. 21 one proba- 
bly would not expect any significant improvement on the mea- 
surement of the neutrino mass over the one obtained using the 
Wmap data alone. This is indeed the case as can be seen from 
Table f2] Our results on Y,m v are in good agr e ement with the 



constraints obtained in Ichikawa et al. (2005); Fukugitaetal 



(2006). Contrary to the results presented in Eisenstein et al. 
(2005J) we do not find any improvement in the measurement 
of Q,(. once the measurement of the low redshift acoustic scale 
is incorporated into the analysis. Although by measuring more 
accurately the value of the Hubble constant one would expect 
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50 60 70 80 90 -1.5 -1 -0.5 0.5 

h q 

Fig. 14. The version of Fig.[T3]with the parameter trio (Ho, O m , 
Weff) replaced by (H , q , jo). 



to better break the geometric degeneracy JBond etal.1 19971) and 
thus should in principle get better constraint on Q.^, this is cur- 
rently not the case. As seen from Fig. [16] by including acoustic 
scale information, model points lie inside the significantly re- 
duced ellipse that is tilted with respect to the bigger one that 
uses the Wmap+HST data alone, in such a way, that the pro- 
jection perpendicular to the flatness line still has practically the 
same width i.e. the constraints on D.^ are also practically iden- 
tical. We can see how the original distribution of the models 
turns towards the and also gets significantly compressed per- 
pendicular to the almost vertical degeneracy line correspond- 
ing to the measurement of the low redshift sound horizon (see 
Eq. (IB. 19b ). In Table [2] along with Y,m v and Q,t we have also 
given a subset of parameters that benefit mostly from the inclu- 
sion of the measurement of the acoustic scale. As expected, the 
parameter constraints are generally getting weaker as we allow 
for more freedom in the models. 



8 
7 
6 

■p 5 

o 


- 

7 




t 4 

3 


' ' ,' S 12 

0.9 




2 
1 


pj^^^^^^^^— 1 


n 




) 0.2 0.4 0.6 0.8 1 

z 



0.2 0.4 0.6 O.i 

z 



1.15 



1.05 




0.95 



0.85 



Fig. 15. Constraints on the low redshift expansion law. Upper 
panel: 2<x credible regions of the look-back time as a function 
of redshift for the Wmap + HST and Wmap + SDSS LRG SHI 
cases. The inset displays these regions after dividing by the 
look-back time corresponding to the best-fit Wmap cosmology. 
Here in addition to the 2cr contours also lcr regions are given. 
Lower panel: Analog of the inset in the upper panel. Here we 
have given the lcr regions as a function of redshift for (starting 
from the bottommost) the Wmap + HST, Wmap + SDSS LRG 
full spectrum, Wmap + SDSS LRG SH2 and Wmap + SDSS 
LRG SHI cases. 



6. Discussion and Conclusions 

In this paper we have performed a MCMC cosmological pa- 
rameter study using the re sults f r om th e recent SDSS LRG 
power spectrum analysis by Hiitsi (2006) along with the CMB 
temperature-temperature and temperature-polar ization angular 



powe r spectra as determin ed by the Wmap team dHinshaw et al 



2003); iKogut et al.l 120031) . We have carried out the analysis 
in two parts: (i) using the Wmap data + the measurement of 
the low redshift sound horizon as found from the SDSS LRG 
redshift-space power spectrum, (ii) using the Wmap data + full 
SDSS LRG power spectrum as shown in Fig. Q] As the for- 
mation of the acoustic features in the large-scale matter distri- 
bution is theoretically very well understood the separate treat- 
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w„„=-1 




Fig. 16. The distribution of models in Q m - Qa plane. Gray 
points correspond to the Wmap+HST case, black ones also in- 
clude the low redshift sound horizon measurement. The dashed 
line represents spatially flat models. Dotted lines provide ap- 
proximate degeneracy directions as given by Eqs. ( IB. 131 ) and 
dB.19t . Dark energy is assumed to be given by the cosmologi- 
cal constant. 



ment for the oscillatory part of the LRG power spectrum is well 
justified. Moreover, in comparison to the full power spectrum, 
which along with the dependence on several cosmological pa- 
rameters requires additional modeling of the redshift-space dis- 
tortions, nonlinear evolution, and biasing L^J the acoustic scale 
depends on only a few cosmological parameters. The CMB 
measurements calibrate the physical scale of the sound hori- 
zon to a very good accuracy. By comparing it with the scale 
inferred from the low redshift LRG power spectrum measure- 
ments, we are able to get a very tight constraint on the Hubble 
parameter: h — 0.708*QlKi if assuming adiabatic initial condi- 
tions, or h — 0.705+QQ37 if additional shift in oscillation phase 
is allowed. Having a tight constraint on h allows us to break 
several parameter degeneracies, and thus helps us to determine 
various parameters like Q,„, Q«/ m , Q.b with a good precision. 
Also, in comparison to the Wmap + HST data, a significantly 
tighter constraint on w e s is obtained. The full results for all 
the parameters are summarized in Table Q] The obtained val- 
ues are in general in a good agreement with several other pa - 
rameter studies e.g. Percival et alJd2002l) : ISpergel et alJd2003l) : 
Tegmark et al.l (120041) . Relatively tight bounds on (Ho, Q m , 
w e ff) or equivalently on (Ho, qo, jo) nei P us to determine the 
low redshift expansion law with significantly higher precision 
than available from the Wmap + HST data alone. If the initial 
fluctuations are constrained to be adiabatic, the measurement of 
the acoustic scale rules out a decelerating Universe, i.e. qo > 0, 
at 5.5<t confidence level. 



In contrast to the acoustic scale measurement, that gave a 
precise value for the Hubble parameter, the full spectrum pro- 
vides us with a good estimate for the shape parameter T = 
D. m h - 0.207 + "['j ^, which is in a very good agreement with 
the one found in [Tegmark et al.l ([2004J). Since in the O m - h 
plane the F = 0,,,/z = const line (see Fig. [T3l is only rela- 
tively weakly tilted with respect to the relevant CMB degener- 
acy direction Q m /i 2 = const, the obtained limits on D. m and h 
are not as strong as the ones obtained from the measurement 
of the acoustic scale. In contrast, the degeneracy lines corre- 
sponding to the low redshift acoustic scale measurement are in 
many cases almost orthogonal to the Wmap + HST "ellipses", 
which explains the stronger constraints for several parameters. 

Throughout most of this work we have focused on spatially 
flat models and also have assumed negligible contribution from 
massive neutrinos. As expected, by relaxing these assumptions 
the bounds on several cosmological parameters loosen signif- 
icantly. The results of this extended analysis for some of the 
parameters are represented in a compact form in Table [2] 

We have also stressed the need to apply cosmological trans- 
formations to the theoretical model spectra before being com- 
pared with the relevant observational spectrum, which is valid 
only in the reference frame of the fiducial cosmological model 
that was used to analyze the data. So far almost all the param- 
eter studies have completely ignored this point, which is prob- 
ably fine for the shallow redshift surveys. On the other hand, 
in case of more deeper surveys like the SDSS LRG, reaching 
z ~ 0.5, these transformations have to be certainly applied. In 
general the line intervals along and perpendicular to the line of 
sight transform differently. Also the transformations depend on 
redshift. We have shown that for the samples like SDSS LRGs, 
with a typical redshift of z ~ 0.35, the single "isotropized" 
transformation taken at the median redshift of the survey pro- 
vides a very accurate approximation to the more complete treat- 
ment. 

For the parameter estimation we have used the SDSS LRG 
power spectrum down to the quasilinear scales, which calls for 
the extra treatment of nonlinear effects, small scale redshift- 
space distortions and biasing. These additional complications 
can be relatively well dealt with t he aid of the Halo Model (see 
AppendixlAl. We have shown in Huts J (120061) that a simple an- 
alytical model with additional four free parameters is able to 
approximate the observed spectrum to a very good precision. 
Also, the Halo Model has been shown to provide a good match 
to the r esults of the sem ianalytical galaxy formation studies 
(see e.g. lCoorav & Shefhll2002l) . In this paper we have shown 
that to a rather tolerable accuracy the above four extra param- 
eter Halo Model spectra (for reasonable values of the parame- 
ters) can be represented as a simple transformation of the lin- 
ear power spectrum with only one extra parameter ( see Fig . [Hi. 



The si milar type of transformation was also used in ICole et al 
d2005h . 



36 Assuming we do not want to exclude the quasilinear scales from 
our analysis. 



In order to investigate the possible biases introduced by the 
specific method used to extract the sound horizon from the 
power spectrum measurements, we have performed a Monte 
Carlo study, the results of what are shown in Fig. [7] 
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Appendix A: Power spectrum from the halo model 

The halo model description of the s patial clustering of galax ies is a 
development of the original idea by iNevman & Scottl Jl952r ). where 
one describes the correlations of the total point set as arising from two 
separate terms: (i) the 1-halo term, that describes the correlations of 
galaxies populating the same halo, (ii) the 2-halo term, which takes 
into account correlations of the galaxies occup ying different halos. 
For a thorough review see ICoorav & Shethl ([2002). Here w e briefly 
give t he results relevant to the current paper (see Seljak 200jJ;|CoorayJ 
l2004h . 

The power spectrum of galaxies in redshift space can be given as: 



P(k) = P l \k) + P 2h (k) , 



(A.l) 



where the 1-halo term: 



P lh (k) 



I dMn 



< m — _, l)m K p {ko-)\u g {k\M)\i' , (A.2) 



V 



1 if (N(N - l)\M) < 1 

2 if (N(N- l)\M)> 1 

and the 2-halo term: 

P 2h (k) = (r> + ^%T g + 1 -tAp\JK) . 

Here: 

<b I i [p\ V^erf(a) 
^ = kcr^-j = —-—, 

T g = I dMn(M)b(M)— ^Ki(kcr)u g (k\M) , 



r v = f 



I dMw 



(MMM^iiko-MMM) . 



(A3) 
(A.4) 

(A.5) 
(A.6) 
(A.7) 
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In the above expressions n(M) is the mass function and b{M) 
halo bias parameter. W e c alculate them usin g the prescription by 
ISheth & TormerJ dl999h and lSheth etal] d200ll) . n represents the mean 
number density of galaxies: 



-/ 



dMn(M)(N\M). 



(A.8) 



We take the mean of the halo occupation distribution in the following 
form: 



(N\M) 



(A.9) 



where Mq and a are free parameters. The second moment is chosen as 
(see Cooray 2004): 



(N(N - 1)|M> = p z (M){N\MY , 

sMil^fe) if M < 10 13 /r>M 



P(M) 



1 



otherwise. 



(A. 10) 



(A.ll) 



/ in Eq. l lA.7b denotes the logarithmic derivative of the linear growth 
factor: / = . ' , u(k\M) and u g (k\M) are the normalized Fourier 
transforms of the dark matter and galaxy density distributions within a 
halo of mass M. In our cal culations we take bot h of these distributions 
given by the NFW profile (Navarro et al. 1997) and the concentration 
parameter c(M) is taken from Bullo ck et al.l ( 120010 . The one dimen- 
sional velocity dispersion of the galaxies inside a halo with mass M is 
taken to follow the scaling of the isothermal sphere model: 



GM 

2-Rvh 



(A. 12) 



where R v i r is the virial radius of the halo and y is a free parameter. 

After specifying the background cosmology the above described 
model has four free parameters: M , a (Eq. l |A.9t ), a (Eq. ( |A. 12| l) and 
M) ow . The last parameter M| ow represents the lower boundary of the 
mass integration i.e. halos with masses below M\ ow are assumed to be 
"dark". 



Appendix B: Fitting formulae for the acoustic 
scales 

The comoving distance traveled by the sound wave since the Big Bang 
up to redshift z can be expressed as: 



*fe) = 



J Cs (z)(l + z')—dz\ 



where the sound speed: 
c 



C S (Z) 



n.z) = 



V3[l+K(z)] 



^,3.04xl0 4 .^! 

4p y Z 



(B.l) 

(B.2) 
(B.3) 



Using t he Friedm ann equation, Eq. l lB.lt can be integrated to yield 
(seee.g. lHi][l995 f"l: 



s(z) = 



In 



3.46 x 10 3 Mpc 

^Cljl 2 ■ Zeq^(Zcq) 



i+ #M 



(B.4) 



37 The result is valid for high enough redshifts as relevant for the 
propagation of the sound waves. 



where the redshift for the matter-radiation equality can be given as: 



2.41 x 10 4 • Q.„,h 2 



(B.5) 



The acoustic scale relevant for the CMB studies is s, = $(z*), where z, 
denotes the recombination redshift. For the "concordance" cosmolog- 
ical model the acoustic scale imprinted in the matter power spectrum 
Sd - s(Zd) is slightly larger than s,. Here Zd denotes the redshift at 
which the baryons are released from the Compton drag of the pho- 
ton field. To find accurate values for z, and Zd one has to carry out a 
full calculation for the recombination history of the Universe. The re- 
sults of these calculations can be conveniently expressed as the follow- 
ing fitting formulas (a ccurate at ~ 1% level) ( IHu & Sugivamal i 1996: 
lEisenstein& Hull 1998b : 

z, = 1048 [l + 0.00124(£Vi 2 r 0738 ] [l + gl (a m h 2 ) S2 ] , (B.6) 

where 

gi = 0.07S3(a b h 2 r 02is [l + 39.5(Ci h h 2 f 163 y' , (B.7) 

g 2 = 0.560 [l+ 21. l(n A /i 2 ) 181 ] _1 , (B.8) 

and 

z d = 1291(fi m /j 2 ) a251 [l+0.659(a m /! 2 ) a828 ] _1 - 

\\ + b,(Q. b h 2 ) b A , 



where 



bi = 0.313(fi„/! 2 r 04W l+0.607(O„ft 2 ) 0674 l 



0.238(fi„^ 2 ) 



(B.9) 

(B.10) 
(B.ll) 



For Q. b h 2 < 0.03 the drag epoch follows the last scattering of the pho- 
tons. 

From the CMB measurements one can determine the angular scale 
9 that corresponds to the sound horizon at decoupling with a good 
accuracy. This angle can be expressed as: 



0: 



s ,h 
dAz.) ' 



(B.12) 



where s, = s(z,), as given in Eq. ( |B.4t , is measured in Mpc, and we 
have added an extra factor of h to the numerator to convert to the 
usual units of h Mpc. Here d ± (z„) is the comoving angular diame- 
ter distance to the last scattering surface, which is strongly dependent 
on the curvature radius R . As s, is only weakly dependent on Q,„h 2 
and Q.bh 2 , it turns out that measurement of 9 is very sensitive to Cl k . 
The dependence of 9 on various cosmological parameters (around the 
"concordance" model point) is given in the upper panel of Fig. IB. II 
Using the following set of model parameters: (Cl„,h 2 , Q.bh 2 , h, C1 DE , 
Weir), the measurement of 9 constrains directly the linear combination 



0.14 I 1 0.73 



■0.45 



h 
0~7l 



or in case of logarithmic variables the combination 
(Ll m h ) (LIde) (h) 



(B.13) 



To measure the sound horizon s d in the large-scale matter distri- 
bution one has to assume some background cosmological model in 
order to convert the observed redshifts to comoving distances. If the 
assumed fiducial model differs from the true cosmology, the measured 
scale will also be distorted. As shown in Chapter 4, for relatively low 
redshift measurements this distortion can be approximated as a single 
transformation for the "isotropized" comoving interval. (In general the 
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Fig. B.l. The dependence of the angular scale corresponding to 
the sound horizon at decoupling (upper panel), and the sound 
horizon as measured from the low redshift matter power spec- 
trum (lower panel), on various cosmological parameters. In the 
lower panel an effective observational redshift z e ff - 0.35 has 
been assumed. The variation of the parameters has been per- 
formed around the "concordance" cosmological model. For 
this model the central values for and s° bs are ~ 0.6° and 

_, d 

~ 107 h Mpc, respectively. 

comoving intervals along and perpendicular to the line of sight trans- 
form differently, and these transformations also depend on redshift.) 
The observed sound horizon s°j s can be approximated as: 



■ hsd 



(B.14) 



where the extra factor of ft is again included to convert to h Mpc, and 
the functions C\\(z), c±(z), which should be evaluated at the effective 
redshift z eff of the observations (e.g. median redshift), are defined as: 

c l|fe) = ^7^ > c ±(z) = -S77T > (B.15) 



E(z) 



<(z) ' 



where 



E{z) 



n r (i + z f + n,„(i + zf + £i*(i + zf + 



n D£ exp I 3[1 + w DE (z')]dz 



dAz) = RoSi 



d H 






d H 





1 sin x if Clk < 


(x)~- 


■Ax if il k = 




( sinh x if il k > , 


C 


r dz' 



(B.16) 



(B.17) 



(B.18) 



The superscript fld refers to the fiducial model. The dependence of sS 
on various cosmological parameters is shown in the lower panel of 
Fig. IB. II The fiducial model here was taken to be the best-fit Wmap 
cosmology and the "true models" were assumed to populate its inter- 
mediate neighborhood. We also assume z eff = 0.35 as in the case of 
the SDSS LRG sample analyzed in this thesis. Then the linear combi- 
nation of parameters constrained by the measurement of sS bs turns out 
to be 



-0.26 



a„/r 



-0.11 



, + 0.94 , 
0.022/ \0.71 



0.14 
or in case of logarithmic variables 



tn J,2-.-0.26 1 - ri j,2\-0.11/j,\0.94/rv \0.12/„ , -.0.15 

(Ll m h ) (Ll b h ) (ft) (Ll DE ) (vf cf f) 



0.71 j \0J3J \-l) 



(B.19) 



Thus, as probably expected, the strongest dependence is on the Hubble 
parameter ft i.e. the precise measurement of the acoustic scale at low 
redshifts should give us a good estimate for ft. 

B.1. The effect of massive neutrinos 

In addition to the dynamical effect that massive neutrinos have on the 
evolution of density perturbations they also lead to the modification of 
the expansion law of the Universe. Assuming three families of mas- 
sive neutrinos with degenerate masses m v one has to replace the term 
n r (l + zf in Eq. iBTTol with: 



0,(1 + z) 4 [1 + 0.1199 • 7(z)] , where : 



Kz) 



f(z) 



j> 



+ f 2 (z)- 



e*+ 1 



dx, 



5967 ■ m v [eV] 

TT~z 



(B.20) 



(B.21) 



(B.22) 



Here the density parameters corresponding to the photons fl y and total 
radiation Q.,. are related by Sl r = 1.68 ■ Q r . 

Since the lower limit for the sum of neutrino masses com- 
ing from the oscillation experiments Em„ > 0.056 eV (e.g. 
Ilxsgourgues & Pastonl20061) is large enough, Eq. JB.21b at low red- 
shifts (as relevant for the acoustic scale measurements using the 
galaxy clustering data) can be very well approximated as: 



Kz) 



1.076 x 10 4 -m v [eV] 

~z 



(B.23) 



i.e. neutrinos act as nonrelativistic matter. In case of the experimen- 
tally lowest allowable in,, this approximation is good to 0. 1 % for red- 
shifts z < 3.3. 

The relative change of the size of the measurable sound hori- 
zon at z = 0.35 as a function of the sum of neutrino masses ~Lm v , 
for the mass range compatible with the constraints arising from the 
Wmap data (Ichikawaet al. 2005; Fukugitaet al. 2006), is shown in 
Fig. IB. 21 The other parameters h ere are ke pt fixed to the best-fit Wmap 
ACDM model values dSpergel et alj2003f) . As we can see the measur- 
able sound horizon is only a weak function of Em r and thus one would 
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Fig. B.2. The relative change of the sound horizon, as measured 
from the low redshift matter power spectrum, as a function of 
the sum of neutrino masses. Here an effective observational 
redshift ZeS - 0.35 has been assumed. The other parameters 
are kept fixed to the best-fit Wmap ACDM model values. 



not expect any strong constraints on neutrino mass from the measure- 
ment of the low redshift acoustic scale. We remind that the relative 
accuracy of the sound horizon measurement as found in iHutsil ( l2006r) 
is in the range 2 - 7%. 



Appendix C: Relation between (Q. m , w e ff) and (<?o, 
70) 

In case of the constant dark energy equation of state parameter w cff , 
deceleration parameter qo and jerk jo at redshift z = can be expressed 
as: 



Qo 



2 2 



n m + 



1 + 2 W <sff(Weff + !) 



Or 



(C.l) 
(C.2) 



For the spatially flat models i.e. fi„, + Q. DE = 1 these relations can be 
uniquely inverted to provide: 



**m 



Weff 



2[;o-go(l +2q )] 
1 + 2(jo - 3go) 

2(3q - jp) - 1 
3(1 - 2q ) ' 



(C.3) 
(C.4) 



